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CN Abstract 

Simplicial Dirac structures as finite analogues of the canonical Stokes-Dirac structure, capturing the topological laws of the 
system, are defined on simplicial manifolds in terms of primal and dual cochains related by the coboundary operators. These 
kj) finite-dimensional Dirac structures offer a framework for the formulation of standard input-output finite-dimensional port- 
i—[ Hamiltonian systems that emulate the behavior of distributed-parameter port-Hamiltonian systems. This paper elaborates on 
the matrix representations of simplicial Dirac structures and the resulting port-Hamiltonian systems on simplicial manifolds. 
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1 Introduction 

A wide class of field theories can be treated as port - 
Hamiltonian systemsP(|van der Schaft fc Maschke, 2002| , 
1 "M. Schoberl fc K. Schlacher^Olip. The Stokes-Dir ac 



structure defined in[]van der Schaft & Maschke, 2002 ) is 
an infinite-dimensional Dirac structure which provides 
a theoretical account that permits the inclusion of vary- 
ing boundary variables in the boundary problem for par- 
tial differential equations. From an interconnection and 
control viewpoint, such a treatment of boundary condi- 
tions is essential for the incorporation of energy exchange 
through the boundary, since in many applications the in- 
terconnection with the environment takes place precisely 
through the boundary. For numerical integration, simu- 
lation and control synthesis, it is of paramount interest 
to have finite-dimesnional approximations that can be 
interconnected to one another. 

Most of the numerical techniques emanating from the 
field of numerical analysis, however, fail to capture the 
intrinsic system structures and properties, such as sym- 
plecticity, conservation of momenta and energy, as well 
as differential gauge symmetry. Mixed finite element 
methods can be constructed in a such a manner that 
a numbe r of important str uctural prope rtie s are pre- 
serve d | |Bossavit, 1998] ), [] ( IHirani, 2003] ), | ( |Hiptmair, 
20021). Most of the efforts have been focused on sys- 



tems on manifolds without boun dary or zero energy 
flow through the boundary. I n [] (Golo, Talasila, van 
der Schaft, & Maschke, 2004) a mixed finite element 



scheme for structure-preserving discretization of port- 
Hamiltonian systems was proposed. The construction is 
clear in a one-dimensional spatial domain, but becomes 
complicated for higher spatial domains. Furthermore, 
the geometric content of the discretized variables re- 
mains moot, in sense that, for instance, the boundary 
variables do not genuinely live on the geometric bound- 
ary. 



Recently in[ (Seslija, van der Schaft, & Scherpen, 2011 1, 
we have suggested a discrete exterior geometry approach 
to structure-preserving discretization of distributed- 
parameter port-Hamiltonian systems. The spatial do- 
main in the continuous theory represented by a finite- 
dimensional smooth manifold is replaced by a homo- 
logical manifold-like simplicial complex and its circum- 
centric dual. The smooth differential forms, in discrete 
setting, are mirrored by cochains on the primal and dual 
complexes, while the discrete exterior derivative is de- 
fined to be the coboundary operator. Discrete analogues 
of the Stokes-Dirac structure are the so-called simplicial 
Dirac structures defined on spaces of primal and dual dis- 
crete differential forms. These finite-dimensional Dirac 
structures offer a natural framework for the formulation 
of finite-dimensional port-Hamiltonian systems that 
emulate their infinite-dimensional counterparts. The 
resulting port-Hamiltonian sys tems are in the standard 
input-output form, unlike in | (Golo, Talasila, van der 
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Schaft, & Maschke, 2004), where the discretized models 
are acausal (given by a set of differential and algebraic 
equations). The explicit input-output form obtained by 
our scheme has the advantage from both numerical and 
c ontrol perspective over the implicit model present ed in 
||( |Golo, Talasila, van der Schaft, fc Maschke, 2004 1. 
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In this paper, we address the issue of matrix repre- 
sentations of simplicial Dirac structures by represent- 
ing cochains by their coefficient vectors. In this man- 
ner, all linear operator from the continuous world can 
be represented by matrices, including the Hodge star, 
the coboundary and the trace operator. Firstly, we recall 
the definition of the Stokes-Dirac structure and port- 
Hamiltonian systems. In the third section, we define 
some essential co ncepts from discrete exterior calculus 
as de ve loped in [] ( Desbrun, Hirani, Leok, fe Marsden^ 
2002 ),[ ( Hirani, 2003[ ) . In order to allow the inclusion of 
non zero boundary conditions on the dual cell complex, 
in[ (Seslija, van der Schaft, & Scherpen, 20111 we have 
adapted a definition of the dual boundary operator that 
leads to a discrete analogue of the integration by parts 
formula, which is a crucial ingredient in establishing sim- 
plicial Dirac structures on a primal simplicial complex 
and its circumcentric dual. We demonstrate how these 
simplicial Dirac structures relate to the spatially dis- 
cretized wave equation on a bounded domain and to the 
telegraph equations on a segment. Towards the end of 
the paper, we consider the existence of structural invari- 
ants, which are crucial for the control by energy shaping. 

Not e. Some preliminary results of this paper are reported 
in[ (Seslija, Scherpen, & van der Schaft, 2012). The ba- 
sic results of Section 3, 4 and 5 are already co ntained 
in[ (Seslija, van der Schaft, & Scherpen, 2011), with a 
main difference that the current presentation does not 
lean onto the heavy nomenclature of algebraic topology, 
but instead emphasizes matrix representations, making 
it more accessible and easier to implement. 

2 Background of port-Hamiltonian systems 



Dirac structures were originally developed in[ (Courant, 
1990] ),|l( |Dorfmaii, 1993b as a generalization of symplec- 
tic, presymplectic and Poisson structures. Later, Dirac 



structures were employed as the geometric formalism un- 
derpinning generalized interconnected and c on strained 
Hamiltonian systems I ( van der Schaft, 2000 ), [ (van der 
Schaft fc Maschke, 2002[ )~ 



2.1 Dirac structures 

Let Q be a manifold and define a pairing on TQ © T* Q 
given by 

«(/i,e 1 ),(/ 2 ,e a ))) = e 1 (/ a ) + e a (/ 1 ). 

For a subspace V of TQ © T* Q, we define the orthogo- 
nal complement V 1 - as the space of all (/i, e 2 ) such that 
(((/i,e 2 ), (/2,e 2 ))) = for all (/ 2 ,e 2 ). A Dirac struc- 
ture is then a subbundle V of TQ © T* Q which satisfies 
V = V^. 



aim is a treatment of open Hamiltonian systems in such 
a way that some of the external variables remain free 
port variables. For that reason, let T b be a linear vector 
space of external flows, with the dual space J-^ of exter- 
nal efforts. We deal with Dirac structures on the product 
space Qx F b . The pairing on (TQ x F b ) © (T* Q x J£) 
is given by 



((/i,/m)> ( e i> e 6,l))' ((/2,/fc, 2 ), (e 2 ,e bj2 ) 
= ei(/ 2 ) + e byl (f b , 2 ) + e 2 (/i) + e bi2 {fb,i)- 



(1) 



A generalized Dirac structure T> is a subbundle of 
{TQ x F b ) © (T*Q x F b ) which is maximally isotropic 
under (1). 

Consider a generalized Dirac structure T> on the prod- 
uct space Q x T b . Let H : Q —> R be a Hamiltonian. 
The port-Hamiltonian system corresponding to a 4- 
tuple (Q, .Ff,, X>, iJ) is defined by a set of smooth time- 
functions {t H> (x{t), f(t),e(t)) e Q x T b x T b \t e I C 
K)} satisfying the equation 

(-x(t),f(t),dH{x{t)),e{t)) eV for tel. (2) 

The equation (2) implies the energy balance ^(x(t)) — 
(dH(x(t)),x(t)) = (e(t)J(t)). 

An important class of finite-dimensional port-Hamiltonian 
systems is given by 

x = J(x)^O) + g(x)e 



f = 9 T (x) 



dH 

dx 



(3) 



where for clarity we have omitted the argument i, and 
J : T* Q — > TQ is a skew-symmetric vector bundle map 
and g : J- b — > TQ is the independent input vector field. 

In this paper we deal exclusively with Dirac structures 
on linear spaces, which can be defined as follows. Let T 
and £ be linear spaces. Given an / G J- and an e G £, 
the pairing will be denoted by (e|/) € R. By symmetriz- 
ing the pairing, we obtain a symmetric bilinear form 
((,)) : Fx£ -> Knaturally given as (((/i,ei), (/ 2 ,e 2 ))) = 
(ei|/ 2 ) + (e 2 |/i>. 

A constant Dirac structure is a linear subspace T> C 
J 7 x £ such that T> = J)- 1 , with _L standing for the or- 
thogonal complement with respect to the bilinear form 



2.2 Stokes-Dirac structure 



The notion of Dirac structures is suitable for the for- 
mulation of closed Hamiltonian systems, however, our 



The Stokes-Dirac structure is an infinite-dimensional 
Dirac structure that provides a foundation for the 
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port-Hamiltonian formulation of a class of distributed 
parameter systems with boundary energy flow I (van der 
Schaft & Maschke, 20021). 



Hereafter, let M be an oriented n-dimensional smooth 
manifold with a smooth (n — 1) -dimensional boundary 
dM endowed with the induced orientation, representing 
the space of spatial variables. By fl k (M), k = 0, 1, . . . , n, 
denote the space of exterior fc-forms on M, and by 
fl k (dM), k = 0, 1, . . . , n — 1, the space of fc-forms on 
dM. 

For any pair p, q of positive integers satisfying p + q = 
n + 1 , define the flow and effort linear spaces by 



n p (M) x n q {M) x n n - p {dM) 



e p q = n n -p(M) x n n - q (M) x n n - q (dM) . 



The bilinear form on the product space F p , q X £ Pjq is 

(({fp ! fq > fb ' e p' e <j' e b)i {fp ! fq : /fc J 6 pJ e g ! e b))) 

p,q £^p,q 

ej A / p 2 + ej A / g 2 + e 2 , A # + e\ A /* 
ej A + eg A 



(4) 



rJJ\/ 



Theorem 1 ([|( |van der Schaft fe M aschk eT2002| )) 

Given linear spaces J- Ptq and£ Ptq , and the bilinear form 
((,)), define the following linear subspace T> of T v ^ q x £ Piq 

^ { {fp: fqi fbi &p: ^q: &b) ^ ^~p.q ^ ^P,<? | 



(5) 



'fp' 
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-(-l) n " 9 tr 
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where d is </ie exterior derivative and tr stands for the 
trace operator on the boundary dM . Then T> = T> , that 
is, T> is a Dirac structure. 

Consider a Hamiltonian density H : fF(M) x fl q (M) — > 
fi"(M) resulting in the Hamiltonian # = / M H E R. 



Setting the flows f p = 



da p r 

dt ' Jq 



3 - and the efforts 



e p = S p H, e q = SqH, where (5 p H,S q H) G a n ~P(M) x 
57™ _9 (M) are the variational derivatives of H at (a p , a q ), 
the distributed-parameter port-Hamiltonian sys- 
tem is defined by the relation 



da v da a 



dt ' 



dt 



,f b ,6 p H,8 q H,e b ) 6 2), f eR. (6) 



3 Basics of discrete exterior calculus 

In the discrete setting, the smooth manifold M is re- 
placed by an oriented manifold-like simplicial complex. 
An n-dimensional simplicial manifold K is a simpli- 
cial triangulation of an n-dimensional polytope \K\ with 
an (n — l)-dimensional boundary. Familiar examples of 
such discrete manifolds are meshes of triangles embed- 
ded in K 3 and tetrahedra obtained by tetrahedrization 
of 3-dimensional manifolds. 



3.1 Cho 



id cochc 



The discrete analogue of a smooth fc-form on the man- 
ifold M is a fc-cochain on the simplicial complex K. A 
fc-chain is a formal sum of fc-simplices of K such that its 
value on a simplex changes sign when the simplex ori- 
entation is reversed. The free Abelian group generated 
by a basis consisting of oriented fc-simplices with real- 
valued coefficients is Ck(K;R). The space Ck(K;R) is 
a vector space with dimension equal to the number of 
fc-simplices in K, which is denoted by N k - The space of 
fc-cochains is the vector space dual of Ck(K; K) denoted 
by C k (K;M) or £1%(K), as a reminder that this is the 
space of discrete k- forms. 



Since ^ = J QM e b A f b , the system is lossless. 



The discrete exterior derivative or the coboundary 
operator d k : VL k d {K) -> n k+1 (K) is defined by du- 
ality to the boundary operator dk+i ■ Ck+i{K;Z) — s- 
Ck(K;Z), with respect to the natural pairing between 
discrete forms and chains. For a discrete form a € (K) 
and a chain Ck+i € Ck+i(K; Z) we define d fc by 

(d fe a,c fc+ i) = (a, (d fe ) T c fe+ i) = (a, d k+1 c k+1 ) , 

where the boundary operator d k +i is the incidence ma- 
trix from the space of (k + l)-simplices to the space of fc- 
simplices and is represented by a spars e N k+ i x N k ma- 
trix containing only or ±1 elements] (Desbrun, Kanso, 
& Tong, 2008). The important property of the bound- 
ary operator is d k o d k +\ = 0. The exterior derivative 
also satisfies d fe+1 o d fc = 0, what is a discrete ana- 
logue of the vector calculus identities curlograd = and 
div o curl = 0. 

3.2 Dual cell complex 

An essential ingredient of discrete exterior calculus is 
the dual complex of a manifold-like simplicial complex. 
Given a simplicial complex K, we define its interior dual 
cell complex -k[K ( block complex i n terminology of al- 
gebraic topology] ( Munkres, 1984[ )) as a graph dual of 
K (gemoetrically) restricted to \K\ . The boundary dual 
cell complex -k^K is a dual to dK. The dual cell complex 
■kK is defined as *K — *iK x -k^K . A dual mesh *\K 
is a dual to K in sense of a graph dual, and the dual of 
the boundary is equal to the boundary of the dual, that 



3 




Fig. 1. A 2-dimensional simplicial complex K and its cir- 
cumcentric dual cell complex *K indicated by dashed lines. 
The boundary of *K is the dual of the boundary of K. 

is d(*K) — -k(dK) — *bK. Because of duality, there is 
a one-to-one correspondence between fc-simplices of K 
and interior (n — fc)-cells of *K. Likewise, to every k- 
simplex on dK there is a uniquely associated (n—l — k)- 
cell on d(*K). Fig. 1 illustrates the duality on a flat 2- 
dimcnsional simplicial complex. 

Two most frequent dualities are the barycentric and cir- 
cumcentric (Voronoi) dualities. An important property 
of the the Voronoi duality is that primal and dual cells 
are orthogonal to each other. This feature dramatically 
simplifies the discrete counterpart of the Hodge star, as 
will be shown in the next subsection. For this reason we 
shall be dealing with the circumcentric duality and re- 
quire that the simplicial complex is well-centered (the 
circumcenters of all simplices of all dimensions lie in the 
interior of the corresponding simplices). 



For more d etails on geometric aspects of these operators 
confer to | (Seslija, van der Schaft, & Scherpen, 20111 
and for an example see Section 6. 

3.3 Discrete wedge and Hodge operator 

There exists a natural pairing, via the so-called primal- 
dual wedge product, between a primal /c-cochain 
and a dual (n — fc)-cochain. Let a k <E QH(K) and 
f3 n ~ k e n™~ k (-kiK). We define the discrete primal- 

n k d (K) x nr k (*iK) -+ 



(a k ,o- k )(P n - k ,*iO- k ) 



dual wedge product A 

R by (a k A f3 n ~ k ,V a k) 
(_ 1 )fc(n-fc)^n-fe A a fe ; y o . fe ) j w here V a k is the n- 
dimensional support volume obtained by taking the 
convex hull of the simplex o~ k and its dual ★id' 1 '. 

The proposed definition of the dual boundary opera- 
tor ensures the validity of the summation by parts rela- 
tion that parallels the integration by parts formula for 
smooth differential for ms. This summation by parts for- 
mula was p resented in[ (Seslija, van der Schaft, & Scher- 
pen, 2011) and its matrix version is given by the next 



proposition. 

Proposition 2 Let K be an oriented well-centered sim- 
plicial complex. Given a primal {k — l)-form a and an 
internal dual (n — k)-discrete form /3; £ fl d ~ k (*{K) and 
a dual boundary form /3b G ^d~ k (*bK), then 



(d^a A A, K) + (-l) k - \a A (dT*A + Jf) 

= (tr^a A(3 h ,dK) . 



Everything that has been said about the primal discrete 
forms carries over to the dual cochains, which can be 
interpreted as covectors. The space of dual fc-cochains 
will be denoted as Vl^ (*iK) . The covectors will be labeled 
by a caret symbol, e.g., (3 £ CljfcK). 

The trace operator tv h : fi*j (K) fl d (dK) is a matrix 
that isolates the members of a fc-cochain vector assumed 
on the geometric boundary dK. 

The dual exterior derivative d™ _fc : Q d ~ k (-k i K) 
n^~ k+1 (*{K) is defined by duality to the primal exterior 
operator d fc as 

d n-k = (_i)fe( d fc-i)T_ 

The negative sign appears as the orientation of the dual 
is induced by the primal orientation. 

The dual boundary exterior derivative d b 1_fc : 



n n d - k (* h K) 



(*iK) is defined as 



2—k 



' = (-l) fc - 1 (tr fe - 1 ) T 



PROOF. For completeness, we give this simple proof. 
Observe that 



(d* -1 aAA,iO = 



(d^c^A 



(-l) fc a T d?-*A 
(-lf{aM^- k f3 h K) 



and 



a A d^ k /3 hl K) = a T d£ A> = ((d^ k ) T a) (3 h 
= (-l) fc - 1 (tr fc - 1 a ) T / 3 b 
= (-l) fc_1 (tr^ 1 a A (3 h ,dK) . 



□ 



The support volumes of a simplex and its dual cell are 
the same, which suggests that there is a natural iden- 
tification between primal fc-cochains and dual (n — k)- 
cochains. In the exterior calculus for smooth manifolds, 
the Hodge star, denoted is an isomorphism between 
the space of Worms and (n — fc)-forms. The discrete 
Hodge star is a map : £l k d (K) — > tt d ~ k (*iK) defined 
by its value over simplices and their duals. In case of 
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the circumcentric duality, the Hodge star *k is a diago- 
nal Nk x Nk matrix with the entry corresponding to a 
simplex a h being |<7 fe |/| *j a k \. 

Another possibility for the construction of the Hodge 
operator is to use Whitney forms. The Whitney map 
is an interpolation scheme for cochains. It maps dis- 
crete forms to square integrable forms that are piecewise 
smooth on each simplex. The Whitney maps are built 
from barycentric coordinate functions and the resulting 
matr ix i s sparse but in g eneral not diagonal [ (Bossavit, 
19981 ), H piptmair, 2002 F - 



The linear operators used in this paper are succinctly 
presented in the following diagram 



Sf>{dK)^—Sf{K): 



: Q. n faK) -<-5 — D. n -\* h K) 



q 1 {8K) < — - — n 1 (K) < = > n™- 1 (*; k) n™- 2 faK) 

tr — 1 



d" 



di 1 



n"- 1 (dK) <^n n -\K) <" n 1 (*iK) -<-^ — n°faK) 

d »-i d P 

v 

Q n (K) < " > Q°(mK) 



((( fpi fq ) fl, e p, e 6 ), (fp, f q , f b , e p , e q , e b )))d 



P,9 



= (el A / 2 + A / 2 + e 2 A /p 1 + e 2 A /* ff) (7) 
+ {elAf* + e 2 b Afi,dK). 



A discrete analogue of the Stokes-Dirac structure is the 
finite-dimen sional Dirac structure constructed in Theo - 
rem 4.1 in[ (Seslija, van der Schaft, & Scherpen, 20111. 
The matrix representation of this Dirac structure is given 
by the following theorem. 



Theorem 3 Given linear spaces J- p q and £p q , and the 
bilinear form ((,)) d - The linear subspace C x £ p 
defined by 



T>d — {(fp, fq, fb 



, 6p, eg, ( 



'fp' 




o (- 


l) r dp 






d n-p 





h 




(-l) p tr™ 


~ p e 1 



V,1 P,Q I 



+ (-l) r 



(8) 



with r = pq + 1, is a Dirac structure with respect to the 
pairing ((,)) d . 

PROOF. Note that since d™~ 9 = (-l)9(d"-P) T and 
d T q = (-l)"~ p (tr"- p ) T , the operator 





d n-p 

(-l)Hr n - p 



(_l)r d n-« (-iyd^ q 






4 Simplicial Dirac structures 

In this section, we develop the matrix representations 
of simplicial Dirac structures. These structures are 
discrete analogues of the Stokes-Dirac structure and as 
such are defined in terms of primal an d duals cochains 
on the underlying discret e manifold [| (Seslija, van der 
Schaft, & Scherpen, 20111. 



The role of the smooth manifold M in the discrete set- 
ting is played by an n-dimensional well-centered oriented 
manifold-like simplicial complex K. The flow and the ef- 
fort spaces will be the spaces of complementary primal 
and dual forms. The elements of these two spaces are 
paired via the discrete primal-dual wedge product. Let 

Ti„ = SI* faK) x (ll(K) x n n d - p {d(K)) 



C d 
-*p,<i 



-n n d - p (K)xni 



- q faK)xn n d - q (d(*K)). 



The primal-dual wedge product ensures a bijective re- 
lation between the primal and dual forms, between the 
flows and efforts. A natural discrete mirror of the bilin- 
ear form (4) is a symmetric pairing on the product space 
Tp q x £p q defined by 



is skew-symmetric, and thus (8) is a Poisson structure 
on the state space Vt p d faK) x fl d (K). □ 

Remark 1 The Dirac structure (8) is purely topological 
and as such does not depend on the choice of the geometric 
duality. Thus, the equivalent result holds in case of the 
barycentric duality. 

The other discrete analogue of the Stokes-Dirac struc- 
ture is defined on the spaces 

fi q = n p d (K) x nifaK) x nr p (d(*K)) 
£ d P .q = nr p faK) x nr q (K) x nr q (dK) . 

A natural discrete mirror of (4) in this case is a symmet- 
ric pairing defined by 

((( fp,fq,fb, &p, e q, e b), (fp , fq , fb , ^p, e q, e b)l) d 



eJ rd 



(Zfd 

fct p,g 



(el A fp + e q A f 2 q + e 2 A f\ + e\ A f\,K) 
+(e 1 b Af b 2 + e 2 b Af b \dK). 
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Theorem 4 The linear space T>d defined by 

T^d, {(/pi fqi fbi e pj e gj e b) G ^p,9 * ^P,q\ 







(-l)P<?+ 1 d™-'>" 


A 




dp* 







+ 














Pb 



(_l)P t r n -%} 



(9) 



is a Dirac structure with respect to the bilinear pairing 

Us- 

PROOF. The simplicial Dirac structure (9) is the 
dual of (8) , and the proof is analogue to that of Theo- 
rem 3. □ 

In the following section, the simplicial Dirac structures 
(8) and (9) will be used as terminus a quo for the geomet- 
ric formulation of spatially discrete port-Hamiltonian 
systems. 

5 Port-Hamiltonian systems 

Let a function H : f^(*i A x ^d( K ) -> R stand for thc 
Hamiltonian (a p ,a q ) H > H(a p ,a q ), with a p € fl^(-k[K) 
and a q G f2 ? (X) . A time derivative of "H along an arbi- 
trary trajectory i 
t G R, is 



(d p (t),a,(t))e^(*ilOx ^(JQ, 



dt n{a * 



dH da z 

A 



9cL 



dt 



dH 

da„ 



da q 
dt 



,K 



(10) 



Remark 2 Introducing a linear negative feedback con- 
trol as e b = (— l)( n ~ p ^ n ~ q * > ~ 1 *b fb, where*h isthe Hodge 
star on the boundary dK, leads to passivization of the 
lossless port-Hamiltonian system, i.e., 4rH < — {fb A 
*t>/&; dK) < 0. Furthermore, if the Hamiltonian is a /Coo 
function with a strict minimum that is a stationary set 
for the system (11), the equilibrium is asymptotically sta- 
ble. A more elaborate control strategy can be the energy 
shaping method as is briefly discussed in Section 8. 

An alternative formulation of a spatially discrete port- 
Hamiltonian system is given in terms of the simplicial 
Dirac structure (9) . We start with thc Hamiltonian func- 
tion (a p ,a q ) H> H(a p ,a q ), where a p G Cl^(K) and 
a q G £l d {-kiK). In a similar manner as in deriving (11), 
we introduce the input-output port-Hamiltonian system 







dt 




da q 




. dt . 









d" 



= (-l) p tr n 





OH 

da p 


+ 







OH 

da q 







, dH 

da„ 



(12) 



In contrast to (11), in the case of the formulation (12), 
the boundary flows fb can be considered to be freely 
chosen, while the boundary efforts e& are determined by 
the dynamics. Note that the free boundary variables are 
always defined on the boundary of the dual cell complex. 



where the caret sign reminds that the quantity lives on 
the dual mesh. The relations between the simplicial- 
Dirac structure (8) and time derivatives of the vari- 
ables are: f p — — -g^S f q — while the efforts are: 



dH 

i)<\ „ 



dU 

da„ 



This allows us to define a time-continuous port- 
Hamiltonian system on a simplicial complex K (and its 
dual *K) by 



d&p 




dt 




da q 




. dt. 




fb 





(-l) r d™ 
d n-p o 

(-l) p tr"- p 





dU 

da p 


+ (-l) r 






dH 







da q 





eb , 



dH 



(11) 



where r = pq + 1. 

The system (11) is evidently in the form (3). It imme- 
diately follows that ^H — (8b A fb, dK), enunciating a 
fundamental property of the system: the increase in the 
energy on the domain \K\ is equal to the power supplied 
to the system through the boundary dK and d(*K). The 
boundary efforts ib are the boundary control input and 
fb are the outputs. 



6 Physical examples 

In this section we consider the discrete wave equation 
on a 2-dimensional simplicial complex and the telegraph 
equations on a segment. 

6.1 Two-dimensional wave equation 

Consider the wave equation fi = —EAu c , with 
u c (t, z) G K, z — (zi,z 2 ) G M, where fi is the mass den- 
sity, E is the Young's modulus, A is the two-dimensional 
Laplace operator, and M is a compact surface with a 
closed boundary. Throughout, the superscript c desig- 
nates the continuous quantities. 

The energy variables are thc 2-dimensional kinetic mo- 
mentum p c , and the 1-form elastic strain e c . The coen- 
ergy variables are the 0-form velocity v c and the 1-form 
stress a c . The energy density of the vibrating membrane 
is H(p, e) = | (e c A a c + p c A v c ), where the coenergy 
and energy variables are related by the constitutive re- 
lations a c — E * e c and v c — l/fi*p c . Thc Hodge oper- 
ator here corresponds to the standard Euclidian metric 
on M. The port-Hamiltonian formulatio n of the vibrat- 
ing membrane in full details is giv en in[ (Golo, Talasila, 



van der Schaft, & Maschke, 2004 



give: 
34). 
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7 V 3 



Fig. 2. A simplicial complex K consists of two triangles. The 
dual edges introduced by subdivision are shown dotted. 

Let us now consider the simplicial Dirac structure under- 
pinning the discretized two-dimensional wave equation. 
The energy variables of the discretized system are cho- 
sen as follows: the kinetic momentum is a dual 2-form 
whose time derivative is set to be / p , the elastic strain 
is a primal 1-form with time derivative corresponding to 
f q , the coenergy variables are a primal 0-form e p and a 
dual 1-form e q . Such a formulation of the discrete wave 
equation is consonant with the simplicial Dirac struc- 
ture (8) for the case when p = n = 2 and q = 1, and is 
given by 



the incidence operator d\ , from the primal edges to the 
primal vertices, is 

[vq-,vx] [v x ,v 2 ] [v 2 ,v ] [vi,v 3 ] [v 3 ,v 2 ] 
v -1 
«i 1 -1 0-1 
v 2 1-10 1 

v 3 1 1 -1 

while the discrete exterior derivative from the vertices to 
the edges is the transpose of the incidence operator, i.e., 
d° = d[. The dual exterior derivative is d 1 = — (d°) , 
while the matrix representation of the d^ operator is 



[V2,Vl\ [Vl,V 3 \ [V 3 ,V4\ [V 4 ,V 2 \ 



"A." 




' -d; 1 " 




6p 






./«. 




d° 












eb , 



tr°e E 



The boundary control variable is the 1-form stress ef,, 
while the output is the boundary velocity. The Hamil- 
tonian of the discrete model is 



1 



1 



. } ^A£*!f+pA - * p,K 



The coenergy variables are the dual 1-form a = ^ 
E *i e and the primal 0-form v = 



*a 1 P- 



The resulting port-Hamiltonian system is 



tip 

di 
di 



d; 1 

-d° 



M*0 





E*! 



1 



fb = -tr *^ 1 ^, 

H 1 



The trace operator is tr = (dj! ) ) T . 

It is trivial to show 

(d% A e q , K) + (e p A (d^e, + d£e 5 ), K) 
= h [h , vi] fb (v ) + e b [«i , t) 3 ] f b (vi ) 
+e b [v 3 ,V4}fb(v 3 ) + e b [v4,v 2 ]fb{v 2 ) , 



(13) 



what confirms that the boundary terms genuinely live 
on the boundary of \K\. 

If the geometric duality is circumcentric, the diagonal 
Hodge operators are 



*o = diag 
*i = diag 



1 



1 



1 



1 



|[*«o]| ' I | ' |[*^2]| ' |[*^3 

I [VQ, Vl}\ \[V1,V2]\ \[v2,Vq]\ [bjV^jj 
|[«l)«o]|' |[U0,*5]|' |[«2,«o]|' |[^3,«5]| 

| [-04, *0 5 ] I 



6.2 Telegraph equations 



where the operators d°, dj, tr° = (d^) 1 ", *i, and * 1 
conform to the diagram at the end of Section 3 when 
n = 2. 



Example. Consider a simplicial complex pictorially 
given by Fig. 2. The primal and dual 2-faces have coun- 
terclockwise orientations. The matrix representation of 



We consider an ideal lossless transmission line on a 1- 
dimensional simplicial complex. The energy variables 
are the charge density q £ Slj(K), and the flux density 

4> £ £lj(*K), hence p = q = 1. The Hamiltonian repre- 
senting the total energy stored in the transmission line 
with distributed capacitance C and distributed induc- 
tance L is 
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'■'0 



V-2 



Vi 



V2n-2 



"2r, 



Li Li L 3 L n L n+ i 

— wl. ^w^—r-JW^'''-rwp^-r-nm^- 

— I— n n — I— rr_ — L 



C 'T *T 



X 



v vi 



l-'3 



L'5 



«2n-l «2ri 



Fig. 3. The primal 1-dimensional simplicaial complex K. By 
construction, the nodes vo and «2n are added to the boundary 
to insure that d(*K) — -k(dK). 



A *n 



.K 



(14) 



where *o and *i are the discrete diagonal Hodge oper- 
ators that relate the appropriate cochains according to 
the following schematic diagram 



Fig. 4. The finite-dimensional approximation of the loss- 
less transmission line when the inputs are voltages and 
the outputs are currents. The inductances L\, . . . , are 
the values that the discrete distributed inductance L takes 
on the simplices [i>o, Vi], ■ ■ ■ , [«2n-i, V2n]', the capacitances 
C*i, . . . C n are the values C takes on [v , v 2 ], . . . , [v 2 

The exterior derivative d° : fl d (K) fl d (K) is the 
transpose of the incidence matrix of the primal mesh. 
The discrete derivative d? : Q, d (*iK ) -> Cl d (*iK) in the 
matrix notation is the incidence matrix of the primal 
mesh. Thus, we have 



Q° d (dK) 4r 



1*0 



1*1 



an 

dq 



— V rep- 



where *b is the identity. 

The co-energy variables are: e p 

resenting voltages and e q 

Selecting f p = — 1| and f q = — ^ leads to the port- 
Hamiltonian formulation of the telegraph equations 





*9- = I currents. 



dq 
' dt 

a± 

' dt 



d° 
d? 



*i A 



fb 



(15) 



tr° *q 1 



where are the input voltages and e;, are the output 
currents. 

In the case we want to have the electrical currents as the 
inputs, the charge and the flux densities would be defined 
on the dual mesh and the primal mesh, respectively. In- 
stead of the port-Hamiltonian system in the form (15), 
the discretized telegraph equations would be in the form 
(11). The charge density is defined on the dual cell com- 
plex as q € Sl d (-kiK) and the discrete flux density is 
4> G tt d (K). The finite-dimensional port-Hamiltonian 
system is of the form 



(d?r 



(d 



0\T 

b) 



-1 1 ••■ 

0-1 1 ••• 

- 
-1 ••• 
- 






-1 1 


1 



(17) 



(18) 



Remark 3 The discr ete analogue of the Stokes-Dirac 
structure obtained m[ (Golo, Talasila, van der Schaft, & 
Maschke, 2004) * s a finite- dimensional Dirac structure, 
but not a Poisson structure. The implication of this on 
the physical realization is that, in contrast to our results, 
the transmission line in the finite- dimensional case is 
not only composed of inductors and capacitors but also 
of transformers. 

The physical realizations of the port-Hamiltonian sys- 
tems (15) and (16) are given on Fig. 4 and Fig. 5, respec- 
tively. Stabilization of either of those systems is easily 
achieved by terminating boundary ports with resistive 
elements, what is a practical application of the passiviza- 
tion explained in Remark 2. 

The accuracy of the proposed method is 1 /n (see [ ( Ses- 
lija, van der Schaft, & Scherpen, 2011)). 

7 Conservation laws 



dq 
' dt 

d± 
' dt 

fb 



d° 
d° 



-1? 
o c 

J L * 1 L J 



,-1 1 



tr u * 



C 



Let us consider the existence of conservation laws and 
eb (16) structural invariants for the port-Hamiltonian systems 

on simplicial complexes. 

7.1 Finite-dimensional invariants 



where 8b are the input currents and fb are the output 
voltages. 



The following proposition gives the conditions for the 
existence of conservation laws in the discrete setting. 
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Fig. 5. The finite-dimensional approximation of the loss- 
less transmission line when the inputs are currents 
and the outputs are voltages. The inductances are: 
Ll = /[«„,.*] L ° = L ([vo,V2]), L 2 = f lV2 ^ L c = L([v 2 ,v4), 
L„ = I {v2n _ 2Mn] L c = L([v2„-2,v 2n ]); the val- 
ues of capacitors are: Ci = Jj- o C c = C([vq,Vi]), 
C > = iHh*] ° C = Cil^M), C Z = C c = C([v 3 ,v 5 ]), 

■ ■ •' C'n+l = /[C 2 „_l,«2n] = ^([ fi 2n-l,W2n])- 

Proposition 5 Consider the port-Hamiltonian system 
(11). Let (a pi a q ) i-> C(d p ,a q ) be a real-valued function. 
Then 



E kerd" 



- — e ker d; 



(19) 
(20) 



jjf C is a conservation law for the port-Hamiltonian sys- 
tem (11) satisfying 



dC 



(21) 



where f£ = -(-l)9(P +1 )tr"- p f£ 



PROOF. Differentiating C along the flow of the system 
(11), we have 



dC 
~dt 



dC da r , dC da n 



da. 



A 



dt da. 



A 



dt 



■ K 



„_.a«V fSc 



— f_-y\q(n-q) ( d n- q 9n\ 

[ ' \ da p 



da„ 



(22) 



dcy on 

da„ 



_(_!)9(P+1) ( tr n-p 



dC 
da„ 



-{-iy q d™ -9 



dC \ dn 

da q J da p 



of n, the result (21) follows iff (19) and (20) hold. □ 

Remark 4 // either e b = or f^ = 0, the quantity 
C satisfying (19) and (20) is a conserved quantity — a 
Casimir function. 

1.2 One- dimensional domain 

An interesting case for which it is possible explicitly to 
solve (19) is when p = n. The matrix d° is nothing but 
the transpose of the incidence matrix d\ , from the set of 
edges to the set of vertices, on a connected graph. It is 
a well-known property of any incidence matrix d\ that 
kerc^ = spanl, where 1 stands for the vector with all 
elements equal 1. A direct consequence of this is that 
1 up to a multiplicative constant. 



§c 



In the one-dimensional case the null space of d? is triv- 
ial, cf. (17), what allows us to explicitly express the con- 
servation law. 

Corollary 6 Consider the port-Hamiltonian system 
(11), with p = q = n = 1, on a one- dimensional 
simplicial manifold given on Figure 3. The quantity 
C p = l T a p = a p ([v ,Vi]) + J2kZi a p ([v 2 k-i,V2k+i]) + 
&p{[v2n-ii v 2n\) satisfies the balance law 



dt 



eb(v ) - e b (v 2n ). 



(23) 



In case of the telegraph equations on the segment M = 
[0, 1], the total charge C q — L q c (t, z)dz as well as the 

total magnetic flux = (j) c (t, z)dz are both con- 
servation laws. In the discrete setting, the only con- 
servation law for the system (16) is the total charge 
C q = 1 T <7 whose derivative along the admissible trajec- 
tories is — efc(vo) — e b (v 2n ). Similarly, the total flux 
C^ = l T (f> in the system (15) satisfies the balance law 

= fb(Vo) - fb(v2n), where f b (v ) and fbihn) are in- 
p ut currents. Thes e result differ from those presented in 
[ (Macchelli, 2011), where both the total flux and total 
charge are conserved. 



8 Energy-Casimir method 

Consider the interconnection of (11) with the (possibly 
nonlinear) integrator 



dC 
dt 

Vc 



dH a 



d( ' 



(24) 
(25) 



(-l)9(d"~P) T 

and d™~ 9 = (— l)"~ p (tr™~ p ) T . Furthermore, regardless 



where we have used the fact that d™ q 



whereC € M ro , g c e R mxNb with N b = dimQ^~ q (d(*K)), 
input u c , output y c , and £ i— > -ff c (C) the controller's 
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Hamiltonian. The interconnection is power-preserving 
with u c = fb and ef, = — y c . The composition is the 
port-Hamiltonian system in the form 



da p 




at 




da q 




at 








. dt _ 







d n-p 








dH cl 
d& p 
dH cl 
da q 
dH cl 

(26) 



with (a p ,a q ,() l ~ >> H c \(a p , a q , (,) is the closed-loop 
Hamiltonian H c \(a p , a q , £) = T-L(a p , a q ) + H C {C,). 

The energy shaping for the system (26) is achieved b y re- 
stricting the beha vior of (26) to a certain subspace] (van 



der Schaft, 2000). To this end, we look at the Casimir 
functions of the closed-loop system. 

Proposition 7 The real-valued junction (a p ,a q ,C) ^ 
C(a p ,a q ,£) is a Casimir function of the closed system 
(26) iff 



dC 
da„ 
dC 

da q 
dC 
9C 



ekerd"- p nker (g c tr n " p ) 



G ker 



d n- q (_l)«- 9d »-« g T 



(27) 



PROOF. Solving ^kC(a p ,a q ,(,) = irrespective of 
H cl directly leads to (27). □ 

Remark 5 Since the structural matrix of the port- 
Hamiltonian system (11) is not of full rank in case when 
g c is identity, not all Casimirs of (26) are of the form 
C(a q ,a q ,() = Si(a p ,a q ) - Q, i = 1, . . . , m. 

Remark 6 In case when p = q = m = n = l and g c = 
[1, 1] the only Casimir for the system (26) is l T a q + (. 



9 Final remaks 

An important avenue for the future research is to make a 
connection between continuous and discretized systems 
in presence of a port-Hamiltonian controller. Here, gaug- 
ing the input-output errors is of a great significance for 
the robustness. 

A possible interconnection of a distributed-parameter 
and a finite-dimensional port-Hamiltonian system can 
be established via the interface supplied by a simplicial 
triangulation of the spatial manifold. Energy shaping of 
the closed-loop system can be understood as a structural 
selection of the finite-dimensional controller. This will 
be explored in the forthcoming publication. 
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